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In this paper we introduce a real-space renormalization transformation for random spin systems 
on 2D lattices. The general method is formulated for random systems and results from merging 
two well known real space renormalization techniques, namely the strong disorder renormalization 
technique (SDRT) and the contractor renormalization (CORE). We analyze the performance of the 
method on the 2D random transverse field Ising model (RTFIM). 



I. INTRODUCTION 



Most physical systems are disordered and the description and modeling of such systems is one of the most special 
problems in condensed matter physics. In the early 70's the role of the disorder in physical systems was discussed 
in several papers. Harris [l| formulated a criterion for the relevance of the weak disorder caused by locally random 
^ ' impurities in the system. According to the criterion, the relevance of the disorder depends on the sign of the critical 
exponent for the specific heat. Just one year later, Imry and Ma formulated another criterion, which points 
out the relevance of the weak disorder in less than four dimensions (the ordered state became unstable against an 
arbitrarily weak random field). As it became clear later these criteria can be understood in terms of the real space 
y I renormalization group (see for instance 3,, JJ and references therein). 

' , ' Real space renormalization group (RG) methods for quantum disordered systems were applied first in the late of 
^ 70's. In the pioneering work on this topic Ma, Dasgupta and Hu considered a spin-1/2 antiferromagnetic Heisenberg 
Q chain, where the coupling strengths were assumed to be stochastically distributed |20|. The authors studied the model 
! at zero temperature using a method that essentially relies on the reduction of the number of degrees of freedom in 
the system. The method attracted a lot of attention, was extensively studied and developed further by Fisher [U 
and has been used to investigate a big variety of systems in one and two dimensions (a review on the real space RG 
^ approach can be found in 9]). The method has recently been referred to as strong disorder renormalization technique 
(N| (SDRT). 

T-H Indeed it turned out that the behavior of a system with randomness is in many cases quite different from the 

non-random case. The main ingredient of the disordered systems, which has no counterpart in the systems without 
' I . disorder, is the existence of so called rare regions, i.e., regions possessing atypical properties (for the phase under 
OO consideration) compared to the rest of the system. It is known that this kind of rare effects can govern the behavior of 
' the systems at long distances and result in exotic phases, e.g. Griffiths-McCoy 0, 0, B [sl phase. It is worth to note 
that the application of the method to two dimensional systems is not straightforward, since it distorts the geometry 
of the underlying lattice, and only numerical calculations are possible. Therefore analytic proofs like asymptotic 
exactness of the SDRT in thermodynamical limit (see [2l|, [l^l) do not apply. We will come back to the SDRT in the 
next section and discuss it in more detail. 

Despite its beauty, the SDRT is a perturbative method and applying it to the finite sized systems may cause a 
■ problem. A non-perturbative RG method, that for our knowledge has not been applied to random systems yet, was 
introduced by Morningstar and Weinstein in [2^. This method is called contractor renormalization (CORE) group 
approach and is especially suited for lattice systems (for CORE applications see [1^, [2^, [13, [H, [2^ ) . By definition, 
this method keeps the eigenvalues of the low energy sector and produces an optimal truncation operator from the 
original Hilbert space to the effective one. In other words, the CORE is a non-perturbative block-spin renormalization, 
which uses exact diagonalization to extract the effective interactions in a coarse grained system. 

Having a non-perturbative method on the one hand and ideas of spatially local renormalization of the system 
from the SDRT on the other hand, we introduce a method that unifies both techniques and is suited to investigate 
two dimensional disordered systems. The purpose of this paper is to show that such merging is possible and results 
in a non-perturbative real space renormalization transformation for 2D quantum systems at zero temperature that 
preserves the underlying lattice geometry. 

A reliable real space renormalization group method for description of the low temperature behavior of some system 
gives information about the long distance properties of the system while keeping the fundamental structure of the 
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ground state. This fact is especially relevant in quantum random systems where the entanglement properties of the 
ground state have been identified as the key feature in understanding the behavior of these materials [H, ■ 

Since the method we are about to introduce involves local real-space renormalization steps we will have to analyze 
errors introduced by these local operations. As a benchmark we are going to use statistical arguments, showing that 
long range interactions are not important in the renormalized system and therefore can be neglected, i.e., if we consider 
a model that initially has only nearest neighbors interactions, we can neglect next nearest or more sophisticated terms 
introduced by renormalization and proceed further with a model that has only nearest neighbor interactions. The 
statistical justification of the fact that our method can be applied locally in the real space and without renormalizing 
the whole lattice at once, is a crucial point of this paper. 

The paper is organized as follows. In the next section we are going to discuss two important real space renor- 
malization techniques in more detail and provide the idea of constructing a novel method for renormalization of 
disordered spins systems on 2D rectangular lattices. In section III. we analyze the performance of the introduced 
renormalization transformation and consider several toy models to prove the negligibility of the long range spin-spin 
interactions that might appear during the renormalization process. We give some open problems and incitations for 
future investigations in outlook section. 



II. REAL SPACE RENORMALIZATION GROUP METHODS AND RANDOM SYSTEMS. 
A. Strong Disorder Renormalization Technique 

The name strong disorder renormalization technique reveals perfectly the idea of real space RG for random systems 
introduced by Ma, Dasgupta and Hu in 20]. 

There are a priori several different situations that can appear in disordered systems in the thermodynamical limit. 
When the size of the system increases and the effective disorder becomes a major effect compared to the thermal or 
to the quantum fiuctuations, this effective disorder can either become 

• smaller and smaller without bound: the system is then controlled by a pure fixed point, 

• larger and larger without bound: the system is then controlled by an infinite disorder or infinite-randomness 
fixed point (IRFP), 

• or it may converge towards a finite level: the system is then controlled by a finite disorder fixed point. 

A class of systems whose critical behavior is governed by an infinite-randomness fixed point (IRFP) is characterized 
by a very broad distribution of couplings and a dynamical exponent z that becomes infinite at the critical point. 
In certain models, any initial disorder, even very small, drives the system towards the IRFP at the large scale: in 
particular^this is the case for the random antiferromagnetic quantum spin chain with local spin systems with S = ^ 
(see also [2l|). 

We will illustrate very briefly a concrete scheme for renormalization of systems with infinitely strong disorder in 
one dimension on example of the random transverse field quantum Ising chain (RTFIC) (for detailed consideration 
see [12]). The system has the following Hamiltonian 

H = -Y,J,,ata]-Y.h,at. (1) 

{ij) » 

The basic strategy is to find the strongest coupling in the chain (it can be either {Jij} or {hi}) and minimize 
the corresponding term in the Hamiltonian. The degrees of freedom associated with this maximum energy scale 
rio = max{ Jij, hi} are then frozen at lower energy scales. 

If the strongest coupling is a field, say hk then the spin ak is put in its local ground state, i.e., in the x-direction, 
causing it to become non-magnetic. Effective interactions are then generated between its nearest neighbors; but, as 
all other nearby couplings are likely to be much smaller than hi, these can be treated by second order perturbation 
theory. This introduces new effective interactions 

J ^ JjkJkj 

•Jij — , , [^) 

where i and j are the nearest neighbors of the k. 

If the strongest coupling is an interaction, say J^;, then two spins are combined forming a cluster which, in the 
zeroth order of perturbation theory, has a double degenerate ground state (both up or both down) and thus can be 
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represented again by an effective two-level particle: a new spin. The effective local magnetic field being applied to 
the cluster (kl) is 

7 hkhi 

hu) - (3) 

which results from the second order perturbation theory, where magnetic fields, acting on two spins are considered to 
be small. 

The magnetization of the cluster will be the sum of magnetizations of single spins k and Z, i.e., it changes additively 
^{ki) = ™/c + ?7^;. Since all new couplings are smaller than the initial one rjp, the energy is rescaled and the maximum 
energy is reduced (for more details see [2l|, H^l). Note that the decimation as described above would change the 
geometry of the system in dimensions higher than D = 1, so that we consider the spins to be the vertices of a 
somewhat random graph with the RG modifying the spatial structure in these larger dimensions [8, §]. 

If the quantum disordered phase is renormalized, the fields eventually tend to dominate the bonds and at small 
values of O almost all decimations are cluster annihilations and the effective interactions connecting them becoming 
weaker and weaker; in the procedure of the renormalization, the system hence becomes a collection of asymptotically 
uncoupled clusters with a broad distribution of effective fields. In the ordered phase, in contrast, the interactions tend 
to dominate the fields at low energies, and most decimations are thus decimations of bonds; eventually this causes an 
infinite cluster to form. The zero temperature quantum transition between these phases is a novel kind of percolation 
with the annihilation and aggregation of clusters competing at all energies at the critical point [1, 0, [ISj El, Ull • 

Before closing this section we would like to point out that the SDRT consists of successive local renormalizations in 
the real space, where no long range interactions are considered. It means that after an elementary renormalization step 
the system is described by a Hamiltonian with only nearest neighbors interactions (if one had started with a nearest 
neighbors interactions Hamiltonian) and no next nearest neighbors appear in the Hamiltonian. Strictly speaking, 
after every RG step the ground state of the effective Hamiltonian will deviate from the ground state of the initial one, 
but the error will become asymptotically small in the thermodynamical limit as has been proven by Fisher in [2^ . 
This is the feature we want to retain in our ansatz later on. 



B. The CORE method 



The CORE is the Hamiltonian version of the Kadanoff- Wilson real space RG transformation for lattice field theories 
and lattice spin systems and relies on contraction and cluster expansion techniques. We briefly sketch the main idea 
of the CORE and how it works and refer to ref. [1^ for details. 

The first step in this method is to choose small clusters, elementary blocks which cover the lattice. After that, 
one picks up some of the clusters (since the CORE was introduced for systems with no disorder and with translation 
symmetry, all clusters are the same) and considers the part of the whole Hamiltonian that corresponds to this cluster. 
In what follows we call this part of the Hamiltonian cluster Hamiltonian. For the cluster Hamiltonian one has to 
choose states that are relevant for the description of physical behavior of the cluster (the number and the form of these 
states can vary depending on the particular model). The span of the chosen states forms the effective Hilbert space of 
the cluster. Then a projection Pg// on the effective Hilbert space of the cluster is constructed. This projection is used 

to obtain the so-called range-1 terms of the Hamiltonian expansion {h^^ = Pef f Hl^iuster^ef f) ■ The range-2 terms 
arise from the Hamiltonian that corresponds to two adjacent (connected) clusters. The states of the effective Hilbert 
space of the connected clusters is obtained by taking tensor products of the states single clusters. Afterwards a unitary 
matrix is constructed by means of which the range-2 terms are produced (this matrix is called triangulation matrix 
[23!|). This procedure is iterated to achieve the range- iV terms. Finally the expansion of the truncated Hamiltonian, 
which is the effective Hamiltonian after single renormalization step is written as 

i <i,j> <i,j,k> 

where /i'^'' stands for range-iV term. For more details and rigorous proof that the truncated Hamiltonian can be 
expanded in this way, we refer the reader to ref. [23j . 

Note that for the construction of the range- terms in the expansion ^ one obtains the eigenvalues {e„} and 
eigenvectors {\n)} by the exact diagonalization of N contiguous clusters. The optimal truncation operator (triangu- 
lation matrix) is obtained by a Gram-Schmidt orthogonalization of the eigenvectors of the Hamiltonian projected on 
the effective Hilbert space. In this way, a basis {|0n)} (the remnant eigenstates of the range- A^ Hamiltonian) is built 
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such that the first vector overlaps with the lowest energy eigenvector and those above, the second one with the second 
lowest and those above and so on, i.e.. 



In fact, this reduced basis stems from the QiJ-decomposition of the overlap matrix between the reduced Hilbert space 
and the space of exact eigenvectors of the complete Hamiltonian [2^ . 

Usually, two situations can occur after several steps of the renormalization. The Hamiltonian either flows to a point 
where it can be solved exactly and the correlation length in the effective lattice model goes to zero, or the system is 
self-similar at every scale, the correlation length diverges and the mass gap goes to zero: at this point, the system is 
said to be at the critical point. 

In summary, the CORE has two major advantages over traditional perturbative real-space renormalization schemes: 

• the CORE is not an expansion in weak/strong bonds between block-spins. Its convergence does not necessarily 
depend on the existence of a large gap to the discarded states of the Hilbert space. 

• the CORE is based on an exact mapping form the original Hamiltonian to an effective Hamiltonian, whose 
truncation error can be estimated numerically by calculating higher orders in the expansion. 

Finishing this section we point out that when the Hilbert space dimension is reduced, the CORE provides a good 
description of the initial states in terms of the renormalized states. In order to estimate the quality of the description 
of the states from the constructed effective Hilbert space, one can use an overlap of the lowest energy states |m) with 
the remnant states when the range-iV term in the expansion is constructed. Note that both issues are related, 

as the closer the truncated space is to the exact one, the smaller is the number of terms that should be kept in the 
cluster expansion for a given error. 



In this section we provide the idea how to construct a real space renormalization method for two dimensional 
disordered systems. The details concerning the accuracy of the method are presented in section [Till 



The real space RG method we are about to introduce combines the SDRT to target the clusters to be decimated 
as the ones with the biggest energy gaps and the CORE as a tool to obtain the effective dynamics at a new scale. 

In the Fig[T]we briefly demonstrate how the renormalization works. The renormalization is applied to a ladder of 
spins which is a) traversing the whole lattice and b) contains the local two spin Hamiltonian with the biggest energy 
gap. Accordingly, the 4x4 square lattice is transformed to become a 4 x 3 square lattice with new effective particles 
and nearest neighbor interactions. The choice to renormalize a whole ladder is made to maintain the lattice geometry 
at every step of the transformation. Indeed the renormalization transformation, which involves a single ladder of the 
initial rectangular lattice, results again in a rectangular lattice and therefore preserves the initial geometry of the 
lattice. Note that the ladder involved in the renormalization transformation is either a column or a row of the lattice. 
In what follows we call the renormalization of the single ladder a renormalization step. 

The choice of the ladder follows the position of the local two spin Hamiltonian with the biggest gap. Once this 
Hamiltonian is found, the whole ladder is renormalized. The criterion of targeting the ladder is arbitrary, but might 
have an impact on the outcome of the procedure for some Hamiltonians. For example one can target the ladder, 
which contains a maximal number of local Hamiltonian with a rather big energy gap, albeit the local Hamiltonians 
with a maximum energy gap does not belong to the ladder. We leave the discussion of the different strategies of the 
ladder targeting as an open question. 

Every renormalization step is a sequence of two basic renormalization transformations. In order to see how these 
basic transformations enter the renormalization, we discuss the renormalization of the ladder in more details. First, 
the ladder is decomposed into four-spin blocks, such that some of the blocks form chains and some of them form 
plaquettes (from FIGlUa) to FIGlUb)). The chain terms correspond to the interaction of every rung (two spins) of 
the ladder to its nearest neighbors. The plaquette terms describe interactions between two rungs inside the ladder. 
Note that every pair of spins in the ladder (the rung) contributes to two plaquette terms and to one chain term. After 
the decomposition, each term, representing one of the two basic lattice substructures is renormalized separately using 
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C. Combining the CORE and the SDRT. 



1. General idea of the method. 
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FIG. 1: (Color online) Renormalization of one ladder from 4x4 lattice. Red dots are new effective 2-level systems, which interact 
via effective couplings between each other (blue dashed lines). Note that the ladder is not necessarily a column, it can be also 
a row of the initial lattice. 
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FIG. 2: (Color online) Four essential constituents of the single renormalization step: a) choice of the relevant ladder b) 
decomposition of the ladder into 4 particle terms: plaquettes and chains c) renormalization of each 4 particle term via CORE 
d) averaging of the local effective terms (dashed circles) and assembling of the effective hamiltonians into a renormalized lattice. 



the CORE. This leads to a set of effective 2 and 3 particle Hamiltonians (FIG. [He)). In the final step the effective 
Hamiltonians are assembled to the renormalized Hamiltonian on the smaller lattice (FIG. [2]d)). 
The renormalization step can hence be summarized as follows: 

1. Target the ladder with the biggest local energy gap. 

2. Define the reduced Hilbert space by the lowest energy sector of every pair of spins in the ladder and the rest of 
the untouched spins in the lattice. 

3. Compute exactly the eigenvalues of the four spin problem (the hardest computational step). 

4. Obtain the Hamiltonian on the next scale and rescale the unit of distance and energy. 

5. Iterate the procedure. 

It is noteworthy that step 2. and 3. of the algorithm rely on an unusual implementation of the CORE method. In 
the introductory part we mentioned that CORE makes use of a uniform blocking of the lattice (elementary blocks 
have the same form because of the translational symmetry and used to construct the range- 1 terms of the Hamiltonian 
expansion). Since we now perform the renormalization transformation locally (translational symmetry does not apply 
in the presence of randomness), we need to introduce a non- uniform blocking. This non- uniform blocking will be 
demonstrated using the example of the renormalization of four spins in a chain configuration. The performance of 
such blocking will be analyzed in the next section. 



2. Elementary steps for the successive renormalization transformation. 



The elementary renormalization transformations of the four spin terms mentioned above can be divided into two 
groups. 
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The first type is renormalization of a plaquette, that results in two new particles and new interaction between them 
(Fig. [3]). We use the CORE to renormalize spins in the plaquette configuration. Each pair forms an elementary 
cluster and used to construct the range-1 term. The interaction between two effective particles is given by the range-2 
term. 




FIG. 3: (Color online) Renormalization of four spins in a plaquette configuration. The renormalization transformation results 
in two new particles 1' and 3' and an interaction between them. 

The second type is a renormalization of 4 spins in a chain configuration (FIG. [3]). The latter introduces effective 
interactions to the neighboring spins that increase the accuracy of the method. To use the CORE as it described 
in FIG. [4] we need to modify it. That is to say the size of the elementary clusters varies. We call this way of 
implementation of the CORE - non-uniform blocking. Two central spins form an elementary cluster as well as 
boundary spins. Constructing the range-1 terms the initial Hilbert space of the boundary spins is kept, whereas the 
effective Hilbert space of the central two spins is spanned by the ground state and by the first excited state of the two 
spin Hamiltonian. The range-2 terms are achieved by constructing the triangulation matrix, while the effective Hilbert 
space is a tensor product of the Hilbert spaces of two boundary spins and the span of the two lowest eigenstates of 
the Hamiltonian of two central spins. This modification reflects the fact that the renormalization transformation is 
local in real space due to intrinsic disorder of the system. We will come back to the non-uniform blocking and analyze 
its performance numerically in the next section. 



1 2 3 4 CORE 1 2' 4 




FIG. 4: (Color online) Renormalization of four spins in a chain configuration. The renormalization transformation results 
in three particles 1, 2' and 4. The interaction between the particles is of the short range character (only nearest neighbors 
interact). 



III. ESTIMATION OF LONG DISTANCE AND MULTI-SPIN INTERACTIONS 

In order to investigate the performance of the elementary renormalization transformations described briefly in 
section Hi C 2l we pick up a particular model that is a 2D random transverse field quantum Ising model (RTFIM). The 
Hamiltonian of the 2D RTFIM possesses Z2-symmetry that can be exploited in the renormalization transformation 
and provides a special form of the effective Hamiltonian after each renormalization step (the same observations were 
made for the ID Ising model in [23]). 

A. Z2-symmetry of the 2D Random Transverse Field Ising Model 

The 2D RTFIM is described by the Hamiltonian 

if = -^J,,a,V|-^/i.af (6) 

(ij) i 

where {Jij} are random interactions and the random transverse fields {hi} leading to the quantum fluctuations. The 
specific form of the distribution will be defined later. 

This Hamiltonian is invariant under the transformation — > —erf (Z2-symmetry). The CORE has to preserve this 
symmetry so that the most general form the renormalized Hamiltonian can take is 

^eff = - E O^''^ = • . . a^j;^ (7) 
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where i is the site index, {/i} = {/ii, . . . , is the multi- index (/Xj G {u, a;, y, ^;}) and the ff{/i}'s are the couphngs. 

Due to the Z2-synimetry of the model the only operators that can appear in the one particle Hamiltonian in the 
cluster expansion are {a'^,(T^}; in the two particle nearest neighbor interactions, the symmetries allow terms of the 
form {a^a^} from the Ising interaction and also {a^a^, ct^ct^} and the only three site operators that can appear are: 
{a^ a'^ , a'' , , a'' , , a'^ , a'^a'^a^ , a^a^a^ , crya^a^ , a^a^a^}. From the above discussion 

we conclude that the Z2-symmetry puts certain constrains on the form of the range- iV terms that can appear in the 
expansion of the effective Hamiltonian 

Exploiting the symmetry arguments we will investigate the relevance of the range-3 and range-4 terms that remain 
in the expansion ([4]), when the renormalization follows the Z2-symmetry 

i <i,j> <i:j,k> 

To achieve the task we will consider several scenarios of non-uniform and uniform blocking in various toy models. 



B. Chain of four spins 

First of all we consider a chain of four spins, which after renormalization becomes a chain of three spins (FIG. [5]). 
(This step is an essential part of renormalization transformation as discussed in section III C 2[) . The encircled pair of 
spins and spins on the boundaries of the chain form the range- 1 Hamiltonians. The effective Hamiltonian consists of 
range-1, -2, and -3 terms. 

• (i i) • • • • 

FIG. 5: (Color online) Renormalization of four spins of the RTFIM in a chain configuration using a non uniform blocking. The 
encircled pair of spins and spins on the boundaries of the chain are used to form range-1 terms for the effective Hamiltonian. 
The circles on the right hand side of the figure correspond to the range-1 terms in the effective Hamiltonian. These circles are 
connected by the lines that correspond to range-2 terms. 

Our goal here is to estimate the range-3 terms, which appear in the effective Hamiltonian. There are 10 possible 




FIG. 6: (Color online) Renormalization of four spins of the RTFIM in a chain configuration using a non-uniform blocking. 
(Left plot) XX coupling between the first and the third particle, {a^ (8> 1 (8> o"^) = 0.041, a{a:c <8) 1 ® (Jx) = 0.704. (Right plot) 
ZZ coupling between the first and the third particle, (a^ 8 1® a^) = 0.003, a{az i® 1 ® ctz) = 0.746. 

terms in the range-3 Hamiltonian that satisfy the Z2-symmetry (see section HH Ap . As our simulations show all this 
terms are negligibly small in the presence of disorder. In FIG. [S] we present the XX (the upper picture) and ZZ (the 
lower picture) couplings between the first and the third particle of the renormalized chain. The initial couplings were 
uniformly distributed on the interval [0, 1] and presented statistics were taken after testing 10^ different configurations. 
As one can see from the FIG.jHlthe resulting distributions of both XX and ZZ couplings are symmetric and centered 
at 0. The standard deviations are 0.704 and 0.746 for XX and ZZ interactions respectively. 
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C. Ladder of six spins. Uniform blocking 



In our next example we consider a ladder of six spins, that we transform to a chain of three spins using the uniform 
blocking (FIG. [7]). The encircled pairs of spins are taken to form range-1 terms in the expansion of the effective 
Hamiltonian. As in the previous example the expansion will comprise up to range-3 terms. 



•7 ^ 



FIG. 7: (Color online) Renormalization of six spins of the RTFIM using a uniform blocking. To form range-1 terms of the 
effective Hamiltonian the encircled pairs of spins are used. The circles on the right hand side of the figure correspond to the 
range-1 terms in the effective Hamiltonian. These circles are connected by the lines, that correspond to range-2 terms. 



In FIG.[5]we present the distributions of XX and ZZ coupling strengths between the first and the third particle in 
the resulting chain. The initial distribution was again a uniform distribution from the interval [0, 1] and we collected 
statistics after testing 10^ configurations. As in the previous example both of the resulting distributions have a peak 
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FIG. 8: (Color online) Renormalization of six spins of the RTFIM using a uniform blocking. (Upper plot) XX coupling between 
the first and the third particle, {ax d?) 1 ® Ux) = 0.064, a{ax (8 1® o-x) = 0.500. (Lower plot) ZZ coupling between the first 
and the third particle, (cr^ (g) 1 ® a^) = 0.0003, a{a^ ® 1 ® a^) = 0.0996. 

by and standard deviations 0.500 and 0.0996 for XX and ZZ interactions respectively. 



D. Ladder of six spins. Non-uniform blocking 

In the last example in this section we consider a ladder of six spins, that one transforms to a plaquette of four spins 
using a non-uniform blocking (see FIG.IH]). Two encircled pairs of spins and two single spins are used to derive the 
range-1 terms of the effective Hamiltonian. In this case the resulting Hamiltonian will contain also range-4 terms. 
Our goal here is to show that range-4 terms present in the effective Hamiltonian can be dropped. 

In FIG. [TO] we present statistics for two of the range-4 terms in the effective Hamiltonian, that satisfy the Z2- 
symmetry of the Ising model. These terms are (TxC^xCxCfx and azcrzo'zcrz- The mean value of both distributions can 
be with a good approximation considered to be zero. The standard deviation is 0.823 and 0.240 for the XXXX and 
ZZZZ term respectively. 

Finally, we present analogous statistics for the corresponding range-3 terms fFIG. [TT|) . These terms correspond to 
the next nearest neighbor interactions in the renormalized model. We can compare these results with the results of 
the previous section, where we considered the transformation of the six spin ladder to a three spin chain. The mean 
value here is 0.034 for XX and 0.032 for ZZ interaction. The corresponding standard deviation is 0.835 and 0.805 
for XX and ZZ interactions respectively. 

From the presented examples we conclude that one can apply the non-uniform blocking to perform the renormal- 
ization transformation locally in the real space. Our numerical results show that the range-3 and range-4 terms are 
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FIG. 9: (Color online) Renormalization of six spins of the RTFIM using a non-uniform blocking. To form range-1 terms of 
the effective Hamiltonian the encircled pairs of spins and two single spins (the not encircled ones) are used. The circles on the 
right hand side of the figure correspond to the range-1 terms in the effective Hamiltonian. These circles are connected by the 
lines, that correspond to range-2 terms. 
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FIG. 10: (Color online) Renormalization of six spins of the RTFIM using a non-uniform blocking. (Left plot) XXXX plaquette 
coupling of the range-4 term, {a^ ^Ux^t^x^iyx) = 0.024, g{X ® X ® X ® X) — 0.823. (Right plot) ZZZZ plaquette coupling 
of the range-4 term, {a^ ® ctz ® fJ^ ® o^) = 0.001, a{ijz ^ (Jz ® cr^ ® a^) — 0.240. 
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FIG. 11: (Color online) Renormalization of six spins of the RTFIM using a non-uniform blocking. (Left plot) XX coupling 
between the first and the third particle, {a^ ®) 1 ® ® 1) — 0.034, a{ax ® 1 ® ® 1) — 0.835. (Right plot) ZZ coupling 
between the first and the third particle, (a^ 8 1 (g) ct^ (g) 1) = 0.032, a{az (g) 1 (g) cr^ 8 1) = 0.805. 



small and average out and therefore can be neglected in further considerations. Indeed, as it can be seen from our 
numerics, there are equal number of couplings with negative and positive signs. These contributions cancel each 
other in average. Since a priori no particular distribution of initial couplings was assumed, this fact substantiates 
the assumption that in the case of appropriate, by means of [23 [. distribution of couplings the contribution from the 
long-range interactions to the effective Hamiltonian become negligible. 

It could, for example, happen that the encircled pair of spins has a non-degenerate ground state and a double 
degenerate first excited state. This is true if for example both local magnetic fields are much stronger than the 
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coupling. Such situation is unfavorable for the construction of a range- 1 term in the expansion of the effective 
Hamiltonian, since in every range-1 term we keep two states. This, as we believe, is the main source of errors that 
cause a rather big variance of the distributions of the strengths of range-3 and range-4 terms presented in this section. 
Now if we assigned a particular coupling strength to each bond and a particular magnetic field to each spin on the 
lattice, we would avoid the error, and the all range-3 and range-4 terms would turn exactly to zero, which is illustrated 
by the fact that all of them have an arbitrarily small mean value. 

In the renormalization transformation, introduced in this paper, one chooses a particular part of a lattice (a ladder), 
that corresponds to a suitable distribution of couplings and magnetic fields. According to previous numerical evidence, 
this choice allows to write the resulting Hamiltonian after every renormalization step in the form, which contains only 
range-1 and range-2 terms 

HoS = -Yl {Jt,^la^ + J5<aJ + Jf^-af aj) - ^ (8) 
with nearest neighbors interactions J and local magnetic field /if. 



E. Renormalization of the basic constituent of the ladder 



In the last part of this section we investigate the performance of the renormalization transformation applied to 
a toy model, which is a basic constituent of a ladder. In other words we investigate the basic constituent of the 
renormalization step, as it is described in section fll C II 




FIG. 12: (Color online) The basic constituent of the renormalization step: two 4-spin chains, whose central spins are coupled. 
The renormalization involves two central spins of both chains and results in two new spin-i particles (red circles) and five 
effective interactions (blue dashed lines). 

The toy model is presented in the FIG. [121 Two chains of four spins are coupled such that the central spins form two 
rungs of a ladder. The renormalization transformation involves two rungs, while the boundary spins keep untouched. 
The basic renormalized systems consists of six particles that interact as shown in FIG. 1121 Red circles correspond to 
effective particles that originate from clustering of two central spins of both chains. 

In the FIG. [13] we compare the spectra of the initial model (with eigenvalues X'^'^'^^) and the model after the 

™ Ij^cxact ^^^^ I 

renormalization (with eigenvalues A„ ) and we define the absolute error as e„ = "^^oxact " ■ The absolute error for 

the first gap is smaller than 10^'^. The error grows slightly, as one considers higher energy levels and is of the order 
of 6 • 10^'^ for the fourth gap. From this observation we conclude that the low energy levels of the initial Hamiltonian 
are reproduced with a very good accuracy. 



IV. CONCLUSIONS & OUTLOOK 



In this paper we have introduced a renormalization transformation for disordered systems on 2D lattices that 
preserves the geometry of the underlying rectangular lattice. The transformation is done using the real space renor- 
malization group method CORE with non-uniform blocking. We tested the ability of the non-uniform blocking on 
the random Ising Hamiltonian. Our numerical tests showed that the physics of the ferromagnetic random Ising model 
can be modeled with nearest neighbor interactions and local magnetic fields in agreement with the conjecture 0]. 

Since Hamiltonian of a spin model can be used to investigate entanglement properties of the model [33| , our method 
provides also a tool for studying the entanglement in 2D disordered quantum spin models. However these investigation 
will be the subject of forthcoming publication. 

Beside that, there are several open problems that can be seen as compendia for future investigations. First, we 
are going to apply the procedure to other observables. We stress again that the presented technique is a cluster 
expansion of the Hamiltonian. Therefore in order to be able to calculate the expectation values of other observables 
{e.g. magnetization) we have to know the corresponding cluster expansion. Second, due to the successive nature 
of the introduced technique (every renormalization step is a sequence of the renormalization transformations, which 
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FIG. 13: (Color online). Renormalization of eight spins of the RTIM using a non-uniform blocking. Histograms of the error in 
the first eigenvalues between the effective hamiltonian and the exact hamiltonian. The mean values of the errors that appears 
in the plots are: first gap (2 10"*), second gap (1.4 10"^), third gap (1.4 10"^), fourth gap (2 10"^) 



involve maximally four spins), there is no exponential growth of the time consumption with a systems size. On that 
account we are going to analyze scaling properties of the introduced method. Third, the way of targeting the region 
that is to be renormalized depends only on the energy spectrum of the local Hamiltonian of two spins. Thus the 
method is applicable to any type of Hamiltonian with nearest neigbbors interactions, whose two lowest lying energy 
states are separated by a gap from the rest of the spectrum, depending on its intrinsic parameters. It is also natural 
to consider spin models with a higher dimensionality. 

Furthermore we argue that there is a rigorous analytical form of the introduced renormalization transformation 
and that the renormalization flow has a certain fixed point. 

Another point of future work is to generalize the method to lattices of non- rectangular geometries. One can consider 
a variety of different lattice structures. The renormalization procedure will depend on the particular form of the lattice 
and must be considered anew in every particular case. 

We close by mentioning that our method offers itself to go beyond the usual randomness and investigate models 
possessing a spin glass phase. 
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